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ABSTRACT 

We analyse the implications of superfluid turbulence for neutron star physics. We 
begin by extending our previous results for the mutual friction force for a straight 
vortex array to account for the self-induced flow which arises when the vortices are 
curved. We then discuss Vinen's phenomenological model for isotropic turbulence, and 
derive the associated (Gorter-Mellink) form for the mutual friction. We compare this 
derivation to a more recent analysis of Schwarz, which sheds light on various involved 
issues. Having discussed isotropic turbulence, we argue that this case is unlikely to be 
relevant for neutron stars. Instead we expect a rotating neutron star to exhibit po- 
larised turbulence, where relative flow drives the turbulence and rotation counteracts 
it. Based on recent results for superfluid Helium, we construct a phenomenological 
model that should have the key features of such a polarised turbulent system. 



1 INTRODUCTION 

Turbulence in superfluid Helium has been exciting experimenters and theorists in low temperature physics for the last 50 
years. Not surprisingly, the literature on the subject is vast, and the research has reached the level where theoretical models 
represent experimental data relatively well (Donnelly 1991). Meanwhile, the possibility that turbulence may be relevant for 
neutron star superfluidity has hardly been discussed at all (although see the recent work by Peralta et al (2005, 2006)). This 
is quite surprising because the standard models for neutron stars are strongly inspired by the two-fluid model for superfluid 
Helium. There may be differences between the two problems, but the underlying multifluid dynamics is essentially the same 
(Mendell 1991a; Prix 2004; Andersson & Comer 2006). One would expect this to be particularly true for the vortex dynamics 
and the processes that lead to vortex tangles and turbulent behaviour. 

The aim of this paper is to discuss turbulence in superfluid neutron stars. In particular, we want to extend the description 
of the mutual friction force (Alpar, Langer & Sauls 1984; Mendell 1991b; Andersson, Sidery & Comer 2006), that couples 
the two components in a superfluid neutron star core, in such a way that we can account for vortex tangles. This is an 
important flrst step along a road that should eventually lead to an improved description of superfluid neutron star dynamics. 
The mutual friction coupling enters in the analysis of three important astrophysics problems. First of all, it sets the timescale 
on which the two core fluids are coupled following a pulsar glitch (or any event that involves a velocity difference) (Alpar & 
Sauls 1988; Jahan Miri 1998; Andersson, Sidery & Comer 2006). Secondly, it is one of the key mechanisms that determine 
the damping rate of neutron star free precession (Jones & Andersson 2001). Finally, the mutual friction is one of the main 
damping agents for superfluid neutron star oscillations. In particular, it has been argued that the mutual friction damping 
completely suppresses the gravitational- wave driven instability of the fundamental f-modes (Lindblom & Mendell 1995). It 
is also one of the key damping mechanisms affecting the analogous instability of the inertial r-modes (Lindblom & Mendell 
2000), see Andersson (2003) for a review. However, all existing results have been obtained by making use of the form of the 
mutual friction flrst suggested by Hall & Vinen (1956) for Helium. Since this model is based on the assumption on a straight 
vortex array it is not relevant when the vortices form a disorganised tangle. The mutual friction force that applies in this case 
is known to be rather different (Vinen 1957). This suggests that, if neutron stars exhibit turbulent behaviour, then much of 
the generally accepted thinking about the interior superfluid dynamics may be wrong. 
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Sx = — where a = ^ ^ , (4) 



2 THE STRAIGHT VORTEX CASE 

We have recently discussed the mutual friction force, and the associated coefficients, that should apply when the vortices are 
straight (Andersson, Sidery & Comer 2006). This should be the relevant case for uniformly rotating systems, provided that 
the vortex array is dynamically stable (which may not bo the case). Our analysis was based on the standard two- fluid model 
(Prix 2004; Andersson & Comer 2006) where the superfluid neutrons, with mass density pn and velocity are distinguished 
from a conglomerate of all charged components, represented by pp and Vp and loosely refered to as the "protons" . In absence 
of dissipation the equations of motion for the two fluids can be written 

+ viVj) [vt + exwf] + Vi ($ + Ax) + Exwf Vit;j = . (1) 
where the contituent indices x and y x can be either n or p. The relative velocity is defined as 

wf'^ = vj - vt . (2) 

Furthermore, 

Mx 1 dE 

Mx = — = — , (3) 
rrix rrix oux 

where E is an energy functional (which represents the equation of state in this formalism) and /ix is the relevant chemical 
potential per unit mass. The coefficients 

2a , dE 

— where a = 7-— 5 

Px OWnp 

encode what is known as the entrainment effect, and f> represents the gravitational potential. For a detailed discussion of 
these equations, see Prix (2004) and Andersson & Comer (2006). Since the equations of motion represent momentum balance, 

we SCO that the momentum per fluid element is 

Pi = m^[vi + e^w^^] , (5) 

This illustrates the main impact of the entrainment. The momentum of the each fluid need not be parallel with that fluids 
transport velocity. This effect is important for neutron star cores because tho strong force endowes each neutron with a cloud 
of protons (and vice versa). This affects the effective mass of the particle and, when it moves, alters the momentum. As 
discussed by, for example, Andersson, Sidery & Comer (2006) the entrainment can be related to the effective proton mass rUp 
via 

2a = ppEp = np{mp — ml) = ripSml . (6) 

Formally, tho above equations represent a macroscopic average over a large number of individual fluid elements (in the 
standard way). This is a natural approach for neutron stars where one is interested in fluid dynamics at lengthscales vastly 
larger than the inter-particle separation. However, in order to assign values to the various parameters (like the entrainment 
a) one must regularly resolve the problem at the microscopic level. This then raises issues concerning the averaging procedure 
which is required to make contact with the macroscopic level. In problems involving vortex dynamics, the situation is made 
even more complex. It is then natural to introduce an intermediate "mesoscopic" level which is sufficiently large that one can 
average over a large collection of particles (and hence discuss "fluid" dynamics in a sensible way) and yet small enough that 
one can resolve the individual vortices. The macroscopic equations then derive from an averaging over a large set of vortices, 
leading to the smooth equations of motion (1). As we will see, this averaging is one of the key problems that must be overcome 
in a description of superfluid turbulence. 

On the mesoscopic level, the motion of a neutron vortex is affected by fluid flow past it (via the Magnus "force"). At the 
same time, the circulation associated with the vortex induces circular flow in the protons because of the entrainment. This 
in turn generates a magnetic field on the vortex, and the scattering of electrons off of this field acts resistively on the vortex 
motion (Alpar, Langer & Sauls 1984). Taking this resistive force to be proportional to the difference in velocity between the 
vortex (wl) and the charged fluid flow we have 

f!=TZ{vf-vl). (7) 

acting per unit length of the vortex. Neglecting the inertia of the vortex core, and balancing /f to the Magnus "force" we 
arrive at the mutual friction acting on the neutrons (per unit length of vortex) 

/mf K>' i k , *j klm ^ 7 np /o^ 

i = B pneijkK-'Wnp+BpnCijke k-'kiWj^ , (8) 

(correcting an overall sign error resulting from Andersson, Sidery & Comer (2006) writing the consitutent indices in the wrong 
order). The same force acts on the vortex and an equal and opposite force affects the protons. The vector k' determines the 
orientation of the vortex (see later) and its magnitude (the quantum of circulation) is k = /i/2mn w 2 x 10~^ cm^/s. Our 
estimated (dimensionless) mutual friction coefficients are 
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B' = B% (10) 

where ccp = Pp/p is the proton fraction, in agreement with previous estimates by Alpar, Langer & Sauls (1984) and Mendell 
(1991b). 

In order to complete the description of the mutual friction force we need to note that the equations of motion (1) are 
relevant on the macroscopic scale. Meanwhile the analysis leading to (8) is mesoscopic, as we are considering a single vortex. 
The macroscopic equations of motion follow if wo average over a large collection of vortices, i.e. wo need to work at a sufficiently 
large lengthscale that such averaging is meaningful. We also need to note that the Magnus effect is already accounted for in 
the left-hand side of (1), which means that we only need to add the smooth averaged form of the force (8) to the right-hand 
side. In the particular case of a straight vortex array, the required averaging is easy. To get the required force per unit volume 
we simply need to multiply (8) by the local number density of vortices per unit area, n„ . 

As we will face this issue again later, it is useful to discuss it in a little bit more detail here. At the fluid dynamics level, 
we introduce the "rotation velocity" in the standard way, 

2ni = e'^'Vjvt , (11) 

This means that, on the macroscopic level we identify (assuming that the lengthscale considered is sufficiently small that we 

can treat Sn as constant) 

e'^'-V.pl = 2m„[nL + e„(np - fi^,)] . (12) 

Meanwhile, at the mesoscopic level the circulation of neutron momentum (not velocity) is quantised. Representing the meso- 
scopic "momentum"^ by p^, we have 

1 ^ _^ h 

rrin 2rnn 

Strictly speaking, we should arrive at the first equation by averaging the second. Denoting the average by angular brackets 
we have (in the case of straight vortices) 

(e'''V,p^) = n.mnK' = e'^'^V^pg (14) 
Hence we see that 

UvK = 2[Qn + en(np - Hn)] . (15) 



ijfcV^Pn = Ki = 7^3— . (13) 



3 ARE SUPERFLUID NEUTRON STARS LIKELY TO BE TURBULENT? 

We want to extend our description of the mutual friction force to the more complicated situation where the vortices are curved 
and may reconnect to form a turbulent tangle. Before we try to do this, it is useful to have some guidance as to whether one 
would expect turbulence to be relevant for superfluid neutron stars. After all, virtually every discussion of neutron star vortex 
dynamics has made the assumption that the vortex array is straight. Should this basic tenet fall, a possibility that has been 
hinted at by, for example, Packard (1972) and Chevalier (1995), we might have to reassess much of our current "wisdom". 

To get an idea of the possible relevance of superfluid turbulence, let us adapt a recent argument due to Finne et al (2003). 
They consider the problem for superfluid '^He, and argue that an intrinsic parameter determines the onset of turbulence. In 
our case, the argument would proceed as follows. First assume that the motion of the "normal fluid" can be neglected, i.e. 
take vf = 0. This is convenient becase we then do not have to worry about the effects of entrainment. In the Helium case it 
may also be quite reasonable because one can assume that the normal fluid is viscous and essentially immobile with respect 
to the walls of the container (at least on a timescale large compared to the viscous timescale). In our case, the equation of 
motion for the neutrons can then be written 

{dt+vi^jK + ^i[...] = -fr' (16) 

Pn 

This is easily rewritten as 



1 We have chosen to refer to the mesoscopic quantity as the "momentum", even though this is not dimensionally correct. As should 
be clear from the discussion we obtain the true momentum by averaging over a unit area. A similar comment applies to the discussion 
following eq. (20). 
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dtv^ + Vi[. . .] = ti^uviuj':, + -}f (17) 

Pn 

where = e'^^V^w^i Taking the curl of this equation, and using (8), we arrive at (assuming that the various parameters are 
constant) 

dttT = (1 - B')V X (tT„ X w„) + Z3V X [w„ x (w„ x v^)] (18) 

Comparing to the standard fluid case, Finne et al (2003) argues that the first term on the right hand side (the inertial term) 
induces the transfer of energy to smaller length scales, i.e. drives the turbulence. Meanwhile, the second term leads to damping 
that serves to stabilise the flow. The relative importance of the two effects is determined by the parameter 

«=T^ (19) 

If q' < 1 one might expect the flow to become turbulent. This is in quite good agreement with the experimental evidence 
provided by Finne et al (2003). They show that turbulence sets in for q < 1.3 or so. Further quantitative evidence for this 
idea has been provided by Henderson & Barenghi (2004). They analyse the spectrum of Kelvin waves in the vortex array, and 
argue that q = 1 represents the crossover point at which these waves change from being overdamped (damping out before 
completing a full oscillation) to propagating. The idea would bo that the vortex oscillations can then lead to reconnections 
with neighbouring vortices. This is obviously more likely to happen if the waves propagate along the vortex array. Anyway, 
in the neutron star case we have already shown that w -C 1. Hence, we would always have g -C 1, and according to 
the above criterion one would expect a superfluid neutron star core to be extremely susceptible to becoming turbulent. Of 
course, there are a number of caveats. Most importantly, the argument does not account for the stabilising effect of rotation. 
Nevertheless, it is clear that we must take the possibility of superfluid turbulence seriously. 



4 ACCOUNTING FOR VORTEX CURVATURE 

The mutual friction force (8) is only relevant when the vortices are straight. This should be the case in a star that is uniformly 
rotating, assuming that the vortex array is dynamically stable. This has generally been taken to be the case for neutron 
stars. However, work on the analogous problem for superfluid Helium indicates that there will exist a critical relative flow 
above which oscillations are induced in the vortices rendering the array unstable (Glaberson, Johnson & Ostermeier 1974). 
This leads to the formation of a complex vortex tangle and a state of turbulence. As a first step towards understanding the 
implications of this for neutron stars we will generalise the mutual friction force to the case of isotropic turbulence. To do 
this, we need to move away from the assumption that the vortices are straight. 

The induced fiow due to the vortex curvature will play a key role in our discussion. The analysis of this problem for 
superfluid Helium (Donnelly 1991) is easily extended to the neutron star case. This is not surprising since the basic nature 
of the rotational vortices is the same. Nevertheless, the derivation is provided in Appendix A. It should be useful for readers 
that are unfamiliar with the Helium results. Moreover, our analysis accounts for the entrainment effect which has (as far as 
wo know) not been considered in previous work. Yet we know that it is important for neutron star dynamics. In fact, one 
would not be at all surprised if (at the end of the day) it were relevant for the Helium problem as well. 

As demonstrated in Appendix A, the vortex curvature induces a contribution to the (mesoscopic) neutron momentum 

?^ = ue,,ks" s'"' (20) 
mn 

where 

and we also know that 

s'i = ki (22) 
Si = k^Vjki (23) 

To get an idea of the value of the parameter i^, which determines the tension in the vortex array, we note that a typical value 
of the size of the vortex core would be ao ~ 100 fm, while we can let b be given by the inter-vortex spacing. This would mean 
that (for the usual triangular vortex configuration) we would have 

b^3.4x 10-' f , , 1 ^ cm (24) 



100 rad/s 
Combining these we see that 



log (A) 

Voo/ 
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2 ° \ 100 rad/s 



In other words, for superfluid rotation rates at the opposite ends of the range of observed pulsar periods, P = 10~^ — 10 s 
(say) the value of only varies by about 30%. It is also worth noting that v is similar in the Helium case. Typical values 
would then be ao ~ 10^* cm and b 10^^ cm, leading to log(6/ao) « 14. 

Having derived the curvature-induced momentum for a single vortex, we can apply the result to the two-fluid model for 
a neutron star core. As discussed by Andersson, Sidery & Comer (2006), the presence of a neutron vortex will affect the flow 
of both neutrons and protons. The induced flow makes a contribution 

^ = vt + e4vf-vT) = ^ (26) 

//tn 1 1 I'll 

§- = vf + e^{vf-vf)=0 (27) 

TTlp 

Solving these equations for the individual velocities, we find 

1 ~ fj ffk /oo\ 

-ftijkS-'s (28) 



1 £n £p 

< = -! vti^ks's'"' (29) 

That is, the difference between the two induced velocities is 

< - < = uti^ks's'"' (30) 
where we have defined 

i>=-—^-—u (31) 

How do these additional contributions to the flow affect the mutual friction force? The answer is quite intuitive. The 

induced velocities simply add to the general flows «x in the expression for the force (8). Adding the velocity difference (30) to 

the overall relative velocity in (8), and noting that 

s x{s xs') = -s" , (32) 
we find that the additional contribution to the mutual friction force can be written 

fr'' = -p^^^i>[B's': + Beijks''s"''] , (33) 
or equivalent ly 

fr'^ = -pnni> [B'k'V.k, + Be.jkk'k'^Vik^] (34) 

To get an expression for the total mutual friction force, these terms should be added to (8). Thus we get the force per unit 
length 

= B'p^eijkK^ w^p + Bpneijke'''"'K^ Kiw'^ - pn!> [B'k^VjKi + BeijkK^ k'"Vik''^ . (35) 

Finally, we determine the macroscopic mutual friction by averaging over a collection of vortices. If we assume that the vortex 

curvature is sufficiently small that there are no vortex intersections, i.e. we ignore possible reconnections and turbulent tangles, 
then the total force is simply arrived at by multiplying the above expression by n„ (as in Section 2). Thus we arrive at the 
final result 

= B' pnUvtijkK' w'lp + BpnTiveijke'''"' k^ Kiw'^ - pr.vriv \B'k^VjKi + BeijkK' k''V ik'']^ . (36) 



5 TURBULENCE 

The curvature contributions to the mutual friction are of key importance for the development of superfluid turbulence. In 

principle, one might expect the above result to remain valid (to some extent) also when the vortices are in a tangle. At least 
as long as the separation between different vortex segments is large enough that the local analysis leading to (say) the induced 
flow remains valid, and as long as it is legitimate to ignore reconnections. However, when the vortices are in a tangle it is 
far from easy to extend the mesoscopic result to the macroscopic fluid dynamics level. The key issue concerns the averaging 
procedure. In the case of a straight vortex array, this issue was trivial. We simply needed to know the number of vortices 
per unit area. Once we allow for a tangle it is no longer obvious how one should average over a collection of vortices. This 
difliculty was recognized already by Feynman in his seminal description of superfluid turbulence (Feynman 1955). One has to 
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be wary of the need to account for vortex loops on a scale smaller than that of the "fluid elements" . These vortex loops need 
to be accounted for even if the "resolution" of the averaged equations is too coarse to identify them. There are a number of 
seriously challenging questions here. 

To make some progress we adopt a phenomenological approach, following the pioneering work of Vinen (1957). We will 
then contrast the results with a more recent analysis due to Schwarz (1982, 1988). The results of the two approaches agree 
qualitatively, even if there are differences at the level of interpretation. 



5.1 Vinen's approach 

In one of his many key papers in this area, Vinen (1957) discussed superfluid turbulence and the associated mutual friction 
force. He argued that, in the case of an isotropic vortex tangle the mutual friction should be proportional to the cube of the 
relative velocity in the flow. This general form for the force had previously been postulated by Gorter & Mellink (1949). In 
this section we will reproduce Vinen's argument for the neutron star problem. This leads to an expression for the mutual 
friction force that can be used when the neutron superfluid is in a homogeneous turbulent state. The analysis is essentially 
phenomenological, but it has the advantage that it leads to an explicit expression for the force in terms of the parameters 
from the straight vortex case, cf. Eq. (8). 

We start from the assumption that the resistive force (7) still applies in the case of a vortex tangle. This seems reasonable 
provided that the vortex segments remain well separated, and we focus on a local part of the vortex^ . Let L be the total vortex 
line length per unit volume and assume that the vortex tangle moves (on average) with the superfluid, i.e. use {v\) « v". In 
the case of isotropic turbulence we then have a force per unit volume (Vinen 1957) 

or 

ff = ^nvf - vf) (37) 

To estimate L, Vinen balances the growth rate due to the Magnus effect (which in turn balances ff in the equation of 
motion for the vortex) to dissipation. To get the growth rate one assumes that it is proportional to (7), i.e. use fe ~ TZwpn = 
PiiKBwpn- Then the fact that it should in addition only depend on pn, L and k, together with a dimensional analysis^ leads to 

dL+ ^ Xlf.L'/' (gg. 

dt PnK 

where Xi is a dimensionless parameter (assumed to be of order unity). Vinen next assumes that the vortex tangle decays in 
the same fashion as ordinary turbulence (i.e. through a Kolmogorov cascade). This leads to an overall decay according to 

dL. ^ (39) 

where X2 is another dimensionless parameter. Combining the growth and the decay, we have 

dL dL+ dL^ XifcL^^'^ X2kL'^ 



dt dt dt PnK 27r 

For obvious reasons, this has become known as "Vinen's equation" . We see that a steady state solution requires 



(40) 



^ ^ (41) 



The final expression for the turbulent mutual friction force is obtained by using this value for L in (37). Hence, we arrive at 



f Sn Pn ( Xi 



2 



/mr ^" l^n I 1 K»3 2 pn /,io\ 

This analysis leads to a mutual friction force that depends on the cube of the relative velocity This behaviour was 
first suggested by Gorter & Mellink (1949), which is why (42) is sometimes refered to as the "Gorter-Mellink" force. It is not 
clear how useful one should expect this representation of the turbulent mutual friction to bo, given the various assumptions 
that went into the analysis. Nevertheless, it has been compared to experimental results in a variety of contexts and appears 
to describe the main phenomena rather well (with typical parameters xi ^ 0.3 and X2 1, see Fig. 2 in Vinen (1957)). This 
indicates that the key assumptions must be "at least approximately true" (Vinen 1957). 



Caveats: It is not clear to what extent this key assumption is justified. One must worry about the role of reconncctions in the vortex 
tangle. It is also not obvious that the resistivity parameters remain as in the straight vortex case. In the neutron star problem one would 
need to consider the local magnetic field induced by entrained protons flowing around the vortex, see Andersson, Sidery &c Comer (2006) , 
and the resultant scattering of electrons flowing through a vortex tangle rather than a straight vortex array. Nevertheless, it is natural 
to make this assumption in this first analysis of the problem. 

^ Note: One could allow additional dimensionless combinations of the various parameters to enter here. Also, because the argument is 
based on a dimensional analysis, it is not easy to see how the dimensionless entrainment would enter. 
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5.2 Schwarz's approach 

In a more recent series of papers, Schwarz (1982, 1988) used numerical simulations to investigate the vortex dynamics. It is 
useful to repeat some of Schwarz's theoretical analysis here, since it highlights the issue of averaging over a vortex tangle. It 
also provides a complement to Vincn's more phcnomcnological discussion. 

Wc start by rccaUing that the complete mutual friction force (per unit length of a vortex) can be written 

/ = KPnB X / X Wnp — US X s"] + KpnB' [s X Wnp — i's"] (43) 

We want to consider an isotropic turbulent state. To do this we focus on a small volume of fluid which contains a vortex 
tangle and through which there is a uniform flow represented by Wnp ■ We will take isotropic to mean that both s* and s" are 
rotationally symmetric with respect to Wnp. We are interested in the mutual friction force per unit volume, V. In principle, 
it follows from the integral 

f^' = ^J fd^ (44) 

where we note that the integrand only has support on the vortex lines. Given the symmetries of the assumed state it is easy 
to see that the terms proportional to B' will not give a net contribution to the integral (Schwarz 1988). We also flnd that 



/ 



S X S X WnprfC — —WnpVLIu (45) 



where we have defined 

h = ^Ji^-4)d^ (46) 
Here we have used Wnp = WnpT^y and 511=/- f\\. To represent the induced flow term we will need 

^'^11 = vhr2 J ^ X ^'d^ (47) 

Combining the two contributions to the force we get 

/"^ = Kp„B[L/||Wp„-i>7iL^''^]f|| (48) 
In these expressions we have used the total line length per unit volume 

^=vjd^ (49) 
Schwarz argues that if L is independent of position and system size then it should scale as 

L = cl{^y (50) 



He further determines the scaling parameter cl by comparing to the numerical simulations. 

Finally, let us simplify the expressions somewhat by assuming that the situation is truly isotropic. Then one would expect* 
/|l =2/3 and /; = 0. This leaves us with the force expression 

j.mf 2 SCl 2 pn 

Ji = l^KPn-^W^pWi (51) 

If we compare this expression to (42) we see that the two expressions would be identical if 



Cl 



(52) 



X2 



This analysis highlights the central role of the averaging, i.e. the integration over all vortex segments in the sample 
volume. In order to be able to analyse less symmetric situations one would have to devise a practical way of carrying out the 
integration. 



* The results of the evolutions carried out by Schwarz (1988) suggest that the isotropic assumption is accurate to 10-20% for /y. At the 
same time Ii 7^ 0. Give the many other uncertainties associated with the neutron star problem, it seems reasonable to assume that the 
turbulent tangle is isotropic here. 



8 Andersson, Sidery & Comer 



6 DIGRESSION: POTENTIAL IMPLICATIONS FOR NEUTRON STARS 

It is interesting to pause at this point and discuss the imphcations that a mutual friction force of form (42) might have for 
neutron star dynamics. It is natural to first compare the relative strength of the alternative mutual friction forces. In the case 
of a straight vortex array we have the Hall-Vinen type force 

/hV » BpnliUvlVpn (53) 

per unit volume. Meanwhile, in the case of isotropic turbulence we have argued that we should use the Gorter-Mellink force 



87r\ (xi 



2 



/gm « ^Pn ^ B'wl^ (54) 



(55) 



3k \^X2 
The relative magnitude of these forces is 

/gm STrVxA's'w'pn 



/hv 3 \^X2, 

Taking Xi/X2 ~ 1 and recalling that n„re « 2Qn, where fin is the rotation rate of the superfluid, we readily arrive at 



/gm _ 47r^g2 _ / r y / Af7\^ / 
/hv ~ 3 Kfin ~ UO«cmy' VOn/ Vlrad/s 



(56) 



rAn = r(flp — fin) (with r the radial distance from the rotation axis) and n — /i/2m„ ~ 
2 X 10~^ cm^/s, together with the typical value S « 4 x 10~*. Finally, as a representative value of the relative rotation rate 
attainable in a neutron star we take An/f2n w 5 x 10~* which is suggested by pulsar glitch data (Lyrie, Shemar & Graham 
Smith 2000). This leads to the final estimate 

(57) 



/hv ~ \W<^cmJ \^5xlO-4y \lsj 



where P is the rotation period. This estimate shows that the relative strength of the two forces depends on the location in the 
star. Near the rotation axis the turbulent mutual friction would be significantly weaker than the force for a straight vortex 
array. Meanwhile, far way from the axis, eg. near the equator in the outer neutron star core, the turbulent force may actually 
be stronger than the standard result. That the relative magnitude of the forces is location dependent is obvious since they 
depend on the relative velocity in different ways. It should be noted that our estimate for the relative strength of the two 
forces differs significantly from that of Peralta et al (2005). They argue that the turbulent force is five orders of magnitude 
weaker than the straight vortex force. The difi^erence is simply due to Peralta et al using a smaller value of the relative velocity 
in their estimate. 

Next let us discuss the implications of fully developed superfiuid turbulence for the dynamical coupling of the two fiuids in 
a neutron star core. We have previously shown (Andersson, Sidery & Comer 2006) that the standard mutual friction coupling 
leads to damping of any relative rotation on a timescale of 

In our discussion, we noted that this estimate is about one order of magnitude smaller than the classic result of Alpar & Sauls 
(1988). We argued that the difference was due to the different scenarios considered. Further consideration of the problem 
shows that, if one imposes global conservation of angular momentum in the Alpar & Sauls (1988) scenario then one arrives 
at a coupling timescale very similar to ours. This point has, in fact, already been made by Jahan Miri (1998). Hence, we 
conclude that our estimate of the coupling timescale is reliable. 

Let us now consider the turbulent mutual friction force. As in Andersson, Sidery & Comer (2006) one can show that the 
relative velocity evolves according to 

ftw°P = (59) 

Xp V^ mg J pn 

In the present case this leads to (with Wnp — > w for clarity) 

dtw w Aw^ — » w = — ^° (60) 
where 

3k \m*J \ x2 
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and Wo = w{t = 0). This behaviour is clearly rather different from that infered in the standard case. In particular, it is clear 
that the coupling is typically much slower than exponential. In fact, it is easy to show that there may be a regime where the 
relaxation is essentially linear (at least locally). This would be the case if 

t < tiinear , where tunear = -r^ (62) 
For typical parameters (the same as in (58)) we find 

,,^.....o-^[^y\^)Q\^y\ ,.3, 

This suggests that a region near the rotation axis of the star (within a metre or so) will relax linearly on a timescale of days. 
This behaviour is certainly sufficiently different from the standard exponential relaxation to warrant further investigation. 



7 POLARISED TURBULENCE 

At this point we have achieved what we originally set out to do. We have generalised the mutual friction force to account for 
the fact that the vortices may not be straight. We have also shown how one can "derive" the force that would be appropriate 
when the vortices form an isotropic tangle and the superfluid exhibits fully developed turbulence. Our results are essentially 
taken over from a number of studies of the analogous problem for Helium. However, as we developed the model it became 
clear that the two cases we have so far considered represent extremes rather than the generic situation. In particular, the 
approach to superfluid turbulence that we adopted is based on "counterflow" experiments in Helium (Donnelly 1991). In this 
context, all vortices are excitations generated by the flow and turbulence sets in above a critical relative velocity. This is in 
sharp contrast to the neutron star problem, where a dense vortex array representing the bulk rotation would be present from 
the outset. In this case the relevance of the turbulent mutual friction force (42) is far from clear. What is clear is that, if one 
assumes fully developed turbulence then there can be no large scale rotation in the superfluid. After all, we assumed that the 
vortex tangle is isotropic. On the macroscopic scale, each "fluid element" would then have to be irrotational. Hence, results 
obtained by simply replacing (8) by (42) are debatable. This casts some doubt on the relevance of the recent simulations 
carried out by Peralta ct al (2005, 2006). On the other hand, it is entirely possible that turbulence will develop locally in a 
neutron star, perhaps in connection with a pulsar glitch. One should certainly consider the implications of this further. Yet 
it is also clear that we need to add another dimension to our analysis. 



7.1 The Helium problem 

In recent work on superfluid Helium, the possible interaction between relative flow and rotation has been considered. It is 

probably fair to say that this problem is still not well understood. There certainly does not yet exist an agreed upon model 
that we can adapt to suit the neutron star problem. Nevertheless, there are useful hints. Much of the recent work was inspired 
by an experiment carried out by Swanson, Barenghi & Donnelly (1983). The experiment concerned counterflow in a rotating 
cylindrical vessel (heated from below) and the results can be summarised as follows. In absence of rotation, turbulence sets 
in above a critical counterflow rate. Representing the relative flow by W, one finds that above a certain critical value Wc the 
vortex line density grows as L ~ W^. Once rotation is considered, it appears that the original Wc is eliminated already at 
very slow rotation rates. Instead, Swanson, Barenghi & Donnelly (1983) find that i) for low values of W there is a flat region 
in L, ii) at a first critical flow Wd there is a slight rise to a second flat region, and iii) beyond a second critical flow rate Wc2 
one sees the anticipated rise with L ~ W^. Nevertheless, rotation appears to add fewer than the anticipated L = Lr = 20,/ k 
lines. 

Simulations (analogous to those of Schwarz (1988)) aimed at understanding the experimental results have been carried 
out by Tsubota, Araki & Barenghi (2003) and Tsubota et al (2004). They confirm much of the phenomenology seen in 
the experiment. The analysis associates the lowest critical velocity Wci with the onset of the so-called Donnelly-Glaberson 
instability (Glaberson, Johnson & Ostermeier 1974). This means that one would have W^i ~ Wdg = 2{2uQ,)^'^. Tsubota et 
al (2004) also observes that, even at high relative flow rates, the disordered tangle has roughly the same "rotation" as the 
initially ordered array. This makes some sense since the overall angular momentum should be conserved. Of course, it also 
suggests that the mutual friction force will not simply be described by (42) in a rotating system. Instead, one should consider 
a polarised turbulent state, where the force is essentially a linear combination of (36) and (42) 

To get some idea of how this may work out, we consider the phenomenological model suggested by Jou & Mongiovi (2004). 
They generalise Vinen's equation by noting that, when both rotation and relative flow are considered then the problem has two 
dimensionless parameters, W/kL}^^ and (Q/kL)^/^. Including terms up to quadratic order, but neglecting in anticipation 
of the relative flow being modest, they write down the following evolution equation for the line density, 
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Figure 1. The experimental results of Swanson, Barenghi & Donnelly (1983), read off from their Figure 3 and represented by crosses 
in the figure, are compared to the phenomenological model of Jou & Mongiovi (2004) discussed in the text. The determination of the 
various parameters in the model is discussed in Appendix B. Clearly, the model (given by the dashed line below Wc2 and the solid curve 
above it) describes the experimental results quite well. 



^ = -/3kL' + {aiW + a2V^)L'^^ 



/3if7 + /32W^\/- L 



(64) 



The coefficients qi and 13 foUow from the original Vinen description, but the other coefficients are new and describe the 
competition between relative flow, which tends to drive the turbulence, and overall rotation, that has a stabilising influence. 
To flnd the stationary state of the system we set dL/dt = 0. This leads to a simple quadratic for L^^^ . This quadratic is 
obviously straightforward to solve in general. Nevertheless, in order to reduce the number of free parameters in the model, 
and simplify the solutions, it is useful to impose the constraint 



fi a\ \ P a-i 
Then the two roots become 



^1/2 



and 



ri/2 ^ ai_W_ , f ^ _ A_ \ _ £^ 



W_^Pi n 

K P2 \ K, 



Analysing the stability of these solutions, one flnds that there exists a critical relative velocity 



Ql 



2^ 

P2 



(65) 



(66) 



(67) 



(68) 



The first solution is stable below this velocity, while the second solution is stable above it. 

This relatively simple model turns out to describe the experimental results quite well. There is a flat region in L below 
Wc and the anticipated quadratic behaviour above Wc- Reading off the experimental data from Figure 3 in Swanson, Barenghi 
& Donnelly (1983) (see also Appendix B) and rescaling to our dimensionless quantities we find that the results for different 
rotation rates are well represented by the model, see Figure 1. This is, in fact, a useful sanity check since it shows that the 
various coefficients do not depend on VF or SI in a significant way. The figure also shows that, taking L 



1/2 



^ — , we have 



L < Lh + Lr. The model cannot account for the small step in L observed at Wc\ but for our present purposes this may not 
be very important. 

Given the apparent success of this relatively simple model it would be interesting to apply it to the neutron star problem. 
Then we first of all need to be able to determine the various new coefficients. To arrive at some reasonable values is relatively 
easy. First compare (64) to (42) to see that 



ai 

P 



2tt 



(69) 
(70) 



Next consider the slow-rotation limit where one would expect the vortex array to be straight. Then L — Lr which means that 
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(32 = V2ai 



(71) 



The final two coefficients are obtained from (65) and the identification of Wc with the critical velocity for the onset of the 
Donnelly-Glaberson instability^ Wdg = 2(2i^a)^/^. This leads to 



/3i = ^-4(-y^' xil3^20 

TT \kJ 



02 = V2 

since jB <C 1. 



(72) 
(73) 



7.2 A phenomenological approach to the neutron star problem 

Now that wc have a basic picture of the combined problem of rotation and relative flow, wc can try to formulate a prescription 
for turbulence in a rotating superfluid neutron star. By necessity, this will require a considerable leap of faith. After all, the 
Helium problem is not completely understood and in that case one has the luxury of being able to carry out laboratory studies. 
What hope do wc have of making progress on the neutron star problem? Wc believe there arc reasons to be optimistic. After all, 
a new perspective on superfluid dynamics is emerging and it will be quite exciting to see if even a very basic phenomenological 
model can help us understand observed phenomena, eg. related to glitches. 

Wc want to work at the macroscopic fluid dynamics level, combining our description of mutual friction for a straight 
vortex array (8) with that for a vortex tangle (42) in a "sensible" way. To do this we will assume that we can average over 
the collection of vortices in each fluid element. In the presence of a turbulent vortex tangle we will assume that the averaging 
can be done on a Icngthscalc large enough that each fluid clement contains an isotropic tangle as well as a "polarised" set of 
vortices that can be represented by a (locally) straight vortex array^. 

Combining what we have learned from the Helium problem, a "workable" strategy might be as follows. First note that 
the turbulent vortex tangle does not contribute to V x Un, while the polarised piece obviously does. Neglecting entrainment, 
which is difficult to account for in this phenomenological model, we let the average 

(V X Vn) = 2ujn = Uvii (74) 

define the vorticity Wn on the macroscopic scale. Taking the absolute value of this we can then identify 

|(V X Vn) \ = 2uJn = Lrk , that is Lr = — |(V X Wn)| , and n„ = Lr (75) 

K 

Wc simply associate the vortex line length due to the local rotation Lr with the density of straight vortices per unit area. 
Finally we make contact with the model discussed in the previous section by taking = uin- 

Next we need the relative flow W in (64). In the Helium experiment this is the counterflow along the rotation ax;is. Hence 
it seems natural to represent it by the projection of the relative flow along the local rotation "axis" . This means that we 

would have 

W = |'u7np ■ l^n\ = |«"np ■ H] (76) 

Here we have assumed that > 0. This is not a physical assumption, in fact, there is no reason why the flow should not be 
in a direction opposite to the rotation axis. However, we need to introduce this restriction to keep contact with (64). 

Now wc have all the ingredients required to work out the total line length L per unit volume. First we should check if 
the relative flow exceeds the critical velocity for the onset of turbulence Wc- If it does not, then the sample is not turbulent 
and it is appropriate to use the Hall-Vinen form for the mutual friction (8). If, on the other hand, W > Wc then we need a 
prescription for working out L. Provided that wo can assume that the system is in "equilibrium" this is easy, and the required 
solution is given by (67). From this solution we readily deduce the total fine length in the turbulent tangle Lt = L — Lr. 
Explicitly, we would have 



/ Ql W\ ai/3i W 1/2 , 
= [ — I H t:;; + 



K 



1 / /3r ' 



4 V /3 



1 



(77) 



Wc arc currently considering this instability in the neutron star case. Our preliminary results suggest that the critical velocity remains 

the same, which would make sense. 

^ Note: This approach may fail unless there is a distinguishable upper limit on the size of the largest vortex loops. Having said that, it is 
clear from the beginning that the separation into straight and tangled vortices in each fluid element is artificial. The assumption is made 
in order to complete the model. To what extent it is sensible depends on future results concerning, for example, post-glitch dynamics. 
At present we are simply trying to develop a model for which "we can do calculations" . 
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This allows us to work out the mutual friction due to the tangle by using Lt in (37). Finally, the contribution from the 
polarised part follows from combining (8) with the deduced vortex density Uv = Lr. The combined mutual friction force for 
the polarised turbulent state (obviously only relevant above the critical velocity) can then be written 

/f f = p^Ln {B'eijkn'w^^ + B£i,fc£'=''"ft^/t,wr - [B'k'Vjm + Beijkn' k'Vik''] } + ^p„/tBwf" . (78) 

This model represents what we think may be the simplest "reasonable" approach to the problem of superfluid neutron 

star turbulence. It has all the key ingredients, combining the standard straight vortex mutual friction with a contribxition from 
an isotropic turbulent tangle. The main additional advantage is that all the steps can actually be worked out in a practical 
situation. Of course, the model is phenomenological and comes with a number of (potentially serious) caveats. Obviously, 
when we write down (77) we assume that the system is in equilibrium. This docs not have to be the case. If one expects the 
timescale for reaching equilibrium to be long, then one can instead assume that L evolves according to (64). The assumption 
that W > may also need to be relaxed. This can be done by considering the directionality in the problem in more detail. 
There have been some efforts to do this in the Helium case, see for instance Nemirovskii & Fiszdon (1995), Lipniacki (2001) 
and Jou & Mongiovi (2005). Common to the ideas discussed in these papers is the use of the Onsager symmetry principle to 
ensure that the second law of thermodynamics is honoured. This is important since the system is dissipative. An interesting 
approach to the problem was outlined in a scries of papers by Gourst (1988, 1989, 1992) and Gcurst & van Boolen (1994, 1995, 
1996). The formalism developed in these papers seems particularly promising given that the equations for the basic two-fluid 
system are derived from a variational principle similar to that used by Prix (2004) . An additional "fluid" , represented by the 
scalar density L and the flow vl, is introduced to describe the turbulent tangle. This seems like a natural approach to take 
and there could be a considerable overlap with the general framework for dissipative multi-fluid systems recently developed 
by two of us (Andersson & Comer 2006). We are planning to investigate this possibility in detail in the future. 



8 CONCLUDING REMARKS 

We have discussed the problem of turbulence and the associated mutual friction force in a superfluid neutron star core. 
Our analysis was heavily influenced by work done on the analogous problem for superfluid Helium over the last 50 years, 
starting with the seminal work of Feynman (1955) and Vinen (1957). Focussing on the aspects of the problem that should 
be relevant for neutron stars, we have derived the form of the mutual friction force that should apply if the system becomes 
completely turbulent. We then argued that the neutron star case is likely to require an understanding of polarised turbulence, 
where rotation has a stabilising influence on the ability of relative flow to induce turbulence. In recognition of this, we 
formulated a phenomenological model that represents a first attempt (for neutron stars) to combine the presence of local 
rotation (represented by a "straight" vortex array") and turbulence (an isotropic vortex tangle). 

These results are intriguing, and possibly quite important. It will certainly be interesting to consider the effects that 
turbulence may have on, for example, pulsar glitch relaxation and neutron star oscillations. Having said that, we appreciate 
that this was our very first attempt at understanding the turbulence problem. The analysis may have taken us some way 
towards the solution to the problem, but more than likely we are only beginning to understand what the real questions are. 



APPENDIX A: THE INDUCED FLOW FOR A CURVED VORTEX 

Let us focus on a single vortex. In other words, we assume that the vortices are well separated, which implicitly puts some 
restrictions on the vortex curvature. To describe the vortex we take Si{^,t), where { is a parameter denoting position along 
the vortex (essentially the arc length), as the position vector of the vortex from an arbitrary origin. The situation is illustrated 
in Figure 2 and it is useful to note that 

/ dsi ^ ,„„, 

s- = ^=«. (79) 

„ _ cPs 

A prime represents a derivative with respect to ^. It is also useful to note the alternative form 

// fi klmr-T I {ci-\\ 

Si = -eijkS-'e ViSm (81) 

To determine the self-induced velocity we must solve for the flow associated with the vortex. To do this we first note 
that it is the circulation of momentum pi that is quantised (as in the main text, we distinguish the mesoscopic momentum 
by a bar), not that of velocity. As discussed in detail by Andersson, Sidery & Comer (2006) this is an important distinction 
to make once we include the entrainment. We need to solve (omitting for the moment any constituent indices) 



2 -3^V,s; (80) 
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Figure 2. A schematic illustration of the variables used to describe the self-induced flow of a superfluid vortex. 

— eijfcV^p* = Ki = (82) 
m 2m 

We can write down a formal solution in terms of a Green's function by introducing a vector potential such that 

ei,fcVM'= = -Pi (83) 
m 

Using the associated gauge freedom to set 

VM,- = (84) 
we see that we need to solve a vector Poisson equation; 

V,-VMi = -Ki (85) 
The formal solution is given by 

M{x) = -1 / -^/y (86) 

47r ] \x - y\ 

From this it follows that 

MS) _ 1 f ^iitEl^d'y (87) 



--/■ 



where W = — y'. To mate use of this solution we note that /t' has support only on the vortex. Thus the volume integral 
collapses to a line integral along the vortex, in terms of ^. Then we have R' = x' — s\ and assuming that the vortex is only 
slightly bent wo can make use of the Taylor expansion 

~ 4 + + (88) 
This also leads to 

ki^s'i + (89) 

Finally, we note that if we want to work out the corrections to the Magnus effect decribed by Andersson, Sidery & Comer 
(2006) we are interested in the induced flow in the vicinity of the vortex. This means that we take the limit — » Soi to arrive 
at 

eijkR'k'' ^~eijks''s"''e (90) 
and 



^^f/^^e (91) 

m 8n J t, 

To carry out the integration we introduce a cut-off at ^ = ao, taken to be the size of the vortex core. This regularises the 
integral and also makes sense since wo have not made any assumptions about the nature of the vortex core. We also introduce 
am upper limit ^ = b, assuming that the contribution to the induced flow can by localised to some part of the vortex. As a 
resonable value for b one can use the intervortex separation. Finally, we arrive at the expression 

Pi ^ ^ f ^ \ fi "k li Ilk /no\ 

= — log — eijfeS-'s =iyeijkS-'s (92) 

m An \aoJ 
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In deriving this result we assumed symmetry with respect to ^ = 0. 



APPENDIX B. DETERMINING THE PARAMETERS IN THE PHENOMENOLOGICAL HELIUM 
MODEL 

The parameters used in Figure 1 were arrived at in the way discussed at the end of section 7.1. Specifically, we first take 
L = Lr below the critical velocity. This fixes 

^ = \/2 (93) 
ai 

A rough fit to the experimental results allows us to determine, cf. (67), 

^ w 5 X 10"^ (94) 
and 

1.23 24.6 (95) 
Wo can now combine these different estimates to get 

^ w 0.071 (96) 
and then 

^ « 1.75 (97) 
Next we assume that the critical velocity Wc2 corresponds to Woo- This means that we should have Wdi = 2^/2^0,, or 
^-2|=2,/^ (98) 
Combining this with the observed critical velocity Wc2 ^ O.USV^ cm/s^^^ we find that 

J « 2.65 (99) 

We now have all parameters in (64) expressed in terms of (3. The values we have deduced agree quite well with those used by 
Jou & Mongiovi (2004). Finally, we learn from Vinen (1957) that X2 ~ 1 which means that we could use 

/3=^«0.16 (100) 

Then all parameters in the model have been specified. 
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